Excitations near the boundary between a metal and a Mott insulator 
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A heterostructure of a semi-infinite metal and a Mott insulator is considered. It is supposed that 
both materials have an identical lattice spacing and hopping integrals and differ in the Hubbard 
repulsion which is negligible in the metal and exceeds the critical value for the Mott transition in the 
insulator. At half-filling and for low temperatures the insulator has the long-range antiferromagnetic 
order. Its low-lying elementary excitations are standing spin waves and a spin-wave mode which is 
localized near the interface and has a two-dimensional dispersion. This mode ejects bulk modes from 
the boundary region. The antiferromagnetic ordering of the insulator induces an antiferromagnetic 
order in the metal where the magnetization decays exponentially with distance from the interface. 
This decay is characterized by the correlation length equal to 5-6 lattice spacings. 

PACS numbers: 73.20.-r, 71.27.+a, 73.40.Ns 



I. INTRODUCTION 

In the last few years an active interest is taken in het- 
erostructures fabricated out of strongly correlated sys- 
tems. Looking for new effects and their possible ap- 
plications a wide variety of systems has been investi- 
gated both experimentally and theoretically. In partic- 
ular, it was established that the interface of Mott and 
band insulators can demonstrate metallic behavior i 1 ' 2 ' 3 ^ 
The similar behavior is expected for the interface of two 
Mott insulators^ With lowering temperature this two- 
dimensional (2D) metal becomes superconducting and 
the transition has some properties of the Berezinsky- 
Kosterlitz-Thouless transition iSil The conduction layer 
can be manipulated by gate voltages which thereby pro- 
vides a way for carrier doping by electric field, in a 
disorder- free way, and is promising for new devices 4 No- 
tice also experiments with ultracold atoms where the co- 
existence of several phases with phase boundaries is often 
observed 

In this paper, we study elementary excitations near the 
boundary of a metal and a Mott insulator which have 
identical simple cubic lattices. The boundary is perpen- 
dicular to one of the crystallographic axes. It is supposed 
that the two crystals have identical hopping integrals and 
differ only in the value of the Hubbard repulsion which 
is vanishingly small in the metal and exceeds the critical 
value U c « 2.8B in the Mott insulator. Here B — 6|t| is 
the halfwidth of the electron band in the considered case 
when only the hopping integral t between nearest neigh- 
bor sites is nonzero. At half-filling the critical value U c 
separates the metallic and the insulating phases in the 
bulk»ii In these conditions, the penetration of the metal- 
lic state into the Mott insulator is negligibly small. 12 
When the temperature is lower than the Neel temper- 
ature the semi-infinite insulator has the long-range anti- 
ferromagnetic order. We investigate its spin excitations 
which are the only low-lying excitations in the consid- 
ered case. The small ratio \t\/U allows us to simplify the 
initial problem, reducing it to the interface between a 
metal and a Heisenberg antiferromagnet. The spectrum 



of the antiferromagnet consists of standing spin waves 
with the dispersion of the bulk three-dimensional (3D) 
spin waves and a mode of surface spin waves which is lo- 
calized within a few layers near the surface and has a 2D 
dispersion. The bandwidth of the surface mode is larger 
than that in the pure 2D case with the same exchange 
constant, though it is smaller than the bandwidth in the 
3D case. The mathematical description of the surface 
mode has much in common with the description of point 
defect states and, as in this latter problem, the localized 
states (the surface mode) eject the bulk states (stand- 
ing waves) from the region near the surface. Then we 
consider the influence of the antiferromagnetic ordering 
of the Mott insulator on the magnetic state of the metal. 
We found that the region of the metal near the boundary 
is also antiferromagnetically ordered. The magnetization 
in the metal decays exponentially with distance from the 
boundary. The correlation length is equal to 5-6 lattice 
spacings and is nearly independent of the ratio \t\/U. 

The description of the considered model and the 
derivation of the effective Hamiltonian are given in 
Sec. II. In Sec. Ill, the eigenvalues and eigenvectors of 
Hamiltonians of the semi-infinite metal and Heisenberg 
antiferromagnet, which are two parts of the effective 
Hamiltonian, are considered. In Sec. IV, the magnetic 
ordering in the metal, which is caused by the antiferro- 
magnetic ordering of the insulator, is studied. Conclud- 
ing remarks are presented in Sec. V. 



II. MODEL 

The axes are chosen in such a way that the metal is lo- 
cated in the half-space l x < —1, while the Mott insulator 
is in the half-space l x > 0. Here l x , l y and l z label sites 
of a 3D cubic lattice and the lattice spacing is set as the 
unit of length. For the half-filled case the Hamiltonian 
of the system reads 



H = t 



E 

<LL'> 



E4. 



2 



u E 



Ux>0 



(1) 



where L = (l x ,ly,lz), 1 = (lyJz), the notation < LL' > 
in the first sum indicates the summation over all pairs of 
nearest neighbor sites, aj^ is the electron creation oper- 
ator with the spin projection a =f or J., n-^ a = a,j jCr ai, cn 
and U is the Hubbard repulsion which, as mentioned 
above, is supposed to be larger than U c and much larger 
than \t\. 

In these conditions the Mott region of the system is an 
insulator with the spin excitations being the only low- 
lying excitations. This fact allows us to simplify essen- 
tially the further consideration by using a unitary trans- 
formation which eliminates terms of the first order in t 
that change the occupation of the Hubbard subbands in 
the Mott region. This unitary transformation is simi- 
lar to that used for the derivation of the t-J Hamilto- 
nian from the Hubbard Hamiltonian. 1 — For this purpose 
it is convenient to switch to the Hubbard operator s 14 > 15 
Xjf = |Li)(Lj| in the Mott region. Here \hi) denotes 
four possible site states - the unoccupied state (i = 0), 
the two singly occupied states (i = a), and the dou- 
bly occupied state (i — 2). The electron creation and 
annihilation operators are connected with the Hubbard 
operators by the relations 

a^=X°r + aX£'' 2 , = X? + aX^~' . (2) 

Here we use an alternative definition of a: a = ±1. In 
these notations Hamiltonian {I) reads 

H = Hq + Hi + H 2 , H a = | W + Xl 2 ) , 

Ux>0 

lacr l m <-l 

+ *E E + a h^ a ^-i^) 

la ; x <-i 

lacr lx>0 

+ { x ?l + i x Z+ x h:^ x uf ( 3 ) 



lcr U>-1 



XcrO vOa | \r2 7 — a v —a,2 \ 



lacr l x >0 



lcr Z,>-1 



where a = (±1,0), (0, ±1) are four unitary vectors con- 
necting nearest neighbor sites in the YZ plane. 

The unitary transformation we are looking for has to 
remove terms of the first order in t which change the 



occupation of the Hubbard subbands in the Mott region 
from the transformed Hamiltonian. Terms of this type 
are collected in the part H 2 in Eq. ([3]). Up to the terms 
of the second order in t the transformed Hamiltonian can 
be written as 

H = e s He~ s 

w H Q + H 1 +H 2 + [S, H ] + [S, H 1 ] + [S, H 2 ] 



+ -[S, [S,H ]]. 
The operator S is looked for in the form 

s = ^EE ff [ x $*,i x C 2 - x hZ x l 

lacr l x >0 



(4) 



ICT l x >0 



A Hx A l,i :c + l ~ A U X A ' 



Ux + 1 



A 10 A L-1 ~~ A 10 A l.-1 



XcrO x^ — cr : 2 x^2, — cr -^Ocr 
1,-1 A 10 ~~ A l,-1 A io 



(5) 



and the parameters £ and £' are determined from the 
condition 



H 2 + [S, H ] = 0, 



(6) 



which eliminates the mentioned terms from the trans- 
formed Hamiltonian (|4]). Using Eqs. ([3]), (O, and (J6]) we 
find 



(7) 



Finally the transformed Hamiltonian reads 

H = *£ 



£ I+a,J„ ( r U« 



lacr (a,<-l 

EE 

lcr J.<-1 



( a l,J.-l )< r°M- 



a U x cr a l,ix-l,cr 



J 



E E S l+a,ix S Hx + J E S l-'x + l S lix 



la !x>0 



l,ix>0 



2JJ2 I SinS 



(8) 



where J = At 2 /U , the components of the spin-i vector 
S L are 5£ = X£ and = i(A^ T -X^ 1 ). In Hamilto- 
nian ijHJ), we have neglected terms of the second order in t, 
which change the occupation of the Hubbard subbands in 
the Mott region or give a correction to the kinetic energy 
in the boundary layer l x = — 1. We neglected also terms 
which describe the intrasubband transport in the Mott 
region. At half-filling such processes are suppressed. 

Notice that the spin bond between the boundary layers 
l x = — 1 and l x — is twice as much the bond in the Mott 
region [in the third term of Hamiltonian ((5J) each bond 
appears twice in the sum]. 
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III. SEMI-INFINITE METAL AND 
ANTIFERROMAGNET 

Let us first consider eigenstates of a semi-infinite metal 
and a Heisenberg antiferromagnet which Hamiltonians 
are contained in Eq. ((SJ. The Hamiltonian of the metal 
reads 

kcr !„<-! 



-4,J B -l,cr a W«<7 + 4i„,<7 a k.**-l,°-) . ( 9 ) 



where the translational invariance of the system in the 
YZ plane was taken into account, 



aid 



e llk ai/ ff , 



is the number of sites in the periodic YZ region, k is 

the 2D wave vector, and — | ^ a elV& - Hamiltonian 
([5]) is diagonalized by the unitary transformation 



fc x 



(10) 



where a/ x fc x satisfies the conditions 

Substituting Eq. (fT0|) into Hamiltonian (|9|) and using con- 
ditions pip we find the following relation for az x fc x : 

a/,+i,fc, + ai x -i,k x = (t~ 1 e kkx - 47 k 2) ) ai xkx (12) 
with the boundary condition 



oti x =o,k x = 0. 



(13) 



Here £kfc x is the eigenvalue of Hamiltonian ([9]). We seek 
the solution of Eq. (fT2")l in the form az x fc x ~ e K ( fex ^ x where 
k has to be real or purely imaginary for the eigenvalue 
of Eq. (fH?)) , 2cosh[«;(fc a ;)], be real. Solutions with real 
k-s do not satisfy condition (|13|) and cannot be mixed 
with solutions with imaginary k-s, since they correspond 
to different energies. Only a linear combination of two 
solutions corresponding to the same eigenenergy (for a 
fixed k) form a new solution. The linear combination of 
two solutions with imaginary and opposite in sign k-s, 
the standing wave, satisfies condition (fT3"|) . 



(14) 



a/xfcx = \/-sin(Mx), 



£kfe x = 2 [cos(k x ) + cos(k y ) + cos(fc z )], 



where k x varies continuously in the range (0, 7r). Thus, 
in Eqs. (|10p and (fTTj) sums over k x have to be substituted 



with integrals and the Kronecker symbol with the Dirac 
delta function. 

Now let us consider the semi-infinite Heisenberg anti- 
ferromagnet described by the Hamiltonian 

H « = ?EE s i+a,i«S Ux + J J2 Sux+iSii.. (15) 



la i x >o 



U«>0 



For low temperatures this system is characterized by the 
long-range antiferromagnetic ordering. Therefore to de- 
scribe its low-lying elementary excitations we use the spin 
wave approximation: 



111 'r-^ 



(16) 



St = P+b 



L °L 



PTb 



L U L' 



S7 = P7 b 



L °L 



where the spin-wave operators &l and b\^ satisfy the Bo- 
son commutation relations and 

n=( 7 r )7 r, 7 r), P± = ~ (l ± e inL ) . 

Substituting Eq. p^|) into Eq. (| 1 5[) and using the trans- 
lation invariance of the Hamiltonian in the YZ plane we 
find 



J E [ 3 ( l - Wo) bljwi 



7 k 
1 



k,Z x >0 
(2) 



- (6ki.&-k,i.+i + fo k 6 -k,ix+i 



(17) 



where 



&k; x = 



iklj 



The usual approach to the diagonalization of bi- 
quadratic forms of the type of Eq. (fT7|) is the use of 
the Bogoliubov-Tyablikov transformation 1 ^ which in the 
present case reads 



b ki x = 



E ( u k; x fc x /?kfc x +fkz x fc x /3L k , fcx ) , 



(18) 



where the operators 0^ k and /3kfc x also satisfy the Bo- 
son commutation relations and therefore the coefficients 
Wki x fc x and fk/ x fc x satisfy the following conditions: 

E ( u k; x fc x u k / ifcx - Vu x k x Vw x k x ) = S l x l' x , 
fe x 

E ( u ki x fc x f-k,Z x fc x - Wk/ x fe x U-k,i x fc x ) = 0, 

A x 

E ( M k'xfe x u ki x fe x - u -k,i x /c x W-k,i x fc x ) 



(19) 



i x >0 



E ("' M ,v; -«k; x fe x u k/xfe J =0. 

; x >0 
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The transformation which is opposite to Eq. (|18|) reads 



«-k,i x fc x 



-k,L 



(20) 



ix>0 



In the new representation Hamiltonian (|17|) is diago- 



nal, 



H a = ^2 E kk jl k Jk kx + const. 
If we use Eqs. fT7 |) -(|50 ]) in the relation 

[Acfc x , H a ] = E^k^ /3kfe x , 

we find equations for the determination of the coefficients 
Ukz x /c x , «kz x fc x and the energy S kfex , 



ElrhU 



kfc x "ki x fc a 



= J 



(2) f 

+27 k V -Ki*kx + r («-k,i x +i,fc x + v-k,i x -i,fcj 



(21) 



-Ekfc.V- 



J 



+ 2 7 k 2) <i xfcx 



3 ( 1- g^ x o ) u-k./,,v, 
1 , . 



2 ,Ul 



t k,i x + l,fe x i "I 

with the boundary conditions 

w k ,i x =-i,fc x = 0> '' M„ i a = 



i,fc x 



(22) 



If, for the time being, we neglect terms proportional to 
<5; x o in Eq. (|2"Tj) . solutions for this set of equations can be 
found in the form 



= ft(fe x )/ x 



where k has to be either real or purely imaginary for the 
energy -E"kfe x be real. Again we find that solutions with 
real k-s do not satisfy boundary conditions (|22|) and can- 
not be admixed to solutions with imaginary k-s, since 
these two groups of solutions corresponds to different 
ranges of energy for a fixed pair (k, — k). Thus, the solu- 
tions of this simplified problem which satisfy boundary 
conditions (|2"2"|) are standing waves 



(23) 



Kl x k x = ^kfe x sm{k x (l x + 1)], 
t'-k,i x fc x = 5 kfex sm[k x (l x + 1)], 
2 3 J + -Ekfe x 



.4 



kfc« 



kfc* 



^ \j('iJ + E^ T ) 2 -(iJ^l)' 



^ J(3J + J B kfcx ) 2 -(3J 7k 3 fc ) x )' 



^kfc x =3J\/1 



(3) 



where 7^ = ^[cos(k x ) + cos(k y ) + cos(fc z )] and k x varies 
continuously in the range (0,7r). Thus, in the above for- 
mulas, summations over k x have to be understood as inte- 
grations in the indicated limits and the Kronecker sym- 
bols of k x have to be substituted with the Dirac delta 
functions. 

The simple exponential solutions are inapplicable, if we 
take into account the previously dropped terms which 
are proportional to <5; x o- To obtain solutions for this 
more complicated problem we use a modification of the 
method applied by I. M. Lifshits for the problem of a local 
defect^ Let us introduce the two-component operator 



B) 



k/ 



?ki x 



and determine the matrix retarded Green's function 



(24) 



where B Mx (t) = e %Hat B Ux e 



iH a t 



with the H a deter- 



mined by Eq. (fl~7|) . Green's function ([24]) satisfies the 
equation 



i-D(ktl x Q = 5(t)6u,f 3 + J 



■> ( 1 - \s lm0 ) r. 



27 k 2) n 



D(ktl x l' x ) + -fi 



D(ktJ x - l,l' x ) + D{kt,l x + l,l' x 



(25) 



where the matrices fi and T3 are given by 



n = 



1 

-1 



T3 = 



1 

-1 



If we define D^(ktl x V x ) as G reen's function which sat- 
isfies Eq. (|25[) without the term proportional to <5; x o, the 
solution of Eq. (125(1 can be written as 

D(ktl x l x ) = b^{Ml x l' x ) 

-- / dt'D^ik,^ ~t,l x 0)D{kt'0l' x ). (26) 
2 J -00 

After the Fourier transformation, 

/oo 
dte iut DQs.tl x l' x ), 
-OO 

we find finally 

D(kul x Q = I)^{kwl x l' x ) - ^D^{kLol x 0) 



T Q + ^b^(kLu00) 



where fp is the 2x2 unit matrix. 



D^(kuOQ, (27) 
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Green's function D^°\\su>l x l' x ) corresponds to the Ha- 
miltonian which eigenvalues and eigenstates are given by 
Eq. ([23jl . This Green's function can be easily calculated, 



D (0 \kujlJ' x ) = / dk x sin[k x (l x + l)]sm[k x (l' x + 1)] 
Jo 

"Acfe. —p, r^Quk 



kfc x 



kfc- 



Qiak^ = 



A 2 



iij 



Ei 



kfe a 



(28) 



ki- r 



where 77 = +0. 

As follows from Eqs. (j27|) and ([28]) , along with the 
standing waves (j2"3"|) Green's function Z^kcj^/^,) may 
have poles corresponding to the surface spin waves. 
These poles are connected with the second term in the 
right-hand side of Eq. (|27|) and they are perceptible in 
the spectral function only near the surface of the an- 
tiferromagnet, since Green's functions D^{ku;l x 0) and 
iM ) (kco0l' x ) in this term decrease rapidly with increas- 
ing the distance l x , l' x from the surface. In Fig. [1] the 
spectral function of spin excitations I~mDn(kujl x l x ) is 
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FIG. 1: The imaginary parts of Green's functions flu (ku)l x l x ) 
(the solid lines) and (k.uil x l x ) (the dashed lines) for l x = 
(a) and l x = 5 (b). k = (0,0. 67r). In part (a), the dash- 
dotted line demonstrates the real part of the denominator in 
the second term in the right-hand side of Eq. (|27|l . 



S 




(0,0) 



FIG. 2: The dispersion of the surface spin- wave mode along 
the symmetry lines of the Brillouin zone (the solid line). 
The dashed line demonstrates the dispersion of the two- 
dimensional spin waves, Eq. (|29|l . 



shown in comparison with the spectral function of the 

standing waves ImD^ (\sll>I x I x ) for different distances 
from the surface. Directly on the surface a pronounced 
peak at ujq w 2.3J and an analogous but less inten- 
sive (for the given k) peak at — a->o are observed in the 
spectrum. From Fig. [IJa) it is seen that the peaks are 
manifestations of the poles of the second term in the 
right-hand side of Eq. (j27|) and therefore correspond to 
the surface mode. Already for l x = 5 the peaks be- 
come indistinguishable in the spectrum and the spectra 
ImDuQsidlxlx) and ImD^ (kwl x l x ) become practically 
identical [see Fig.QJb)]. 

The similarity in the description of impurity states 
and the surface mode leads to some similarity in their 
properties. As seen from Fig. [2a), nearly all quasipar- 
ticle weight of the spectrum is concentrated in the peak 
of the surface mode, while bulk modes manifest them- 
selves as a weak shoulder to its right - as in the impurity 
problem,— localized states (the surface mode) eject bulk 
states (standing waves) from the region near the surface. 

Figure [2] demonstrates the dispersion of the sur- 
face spin waves derived from the spectral function 
lmDn(kujl x l x ) for l x = 0. The dispersion (the solid line) 
is shown along the symmetry lines of the Brillouin zone. 
The shape of the curve resembles the dispersion of the 
2D spin waves^ 



w k = 2J\ 1 



(2) 

7 k 



(29) 



shown with the dashed line in the figure. However, 
the frequencies of the surface excitations are somewhat 
higher, though their bandwidth is smaller than the band- 
width 3 J of the 3D spin waves [see Eq. (|23l) ]. 

In the zeroth order values of the spin operators S^-i 
on the metal side of the considered heterostructure van- 
ish. Therefore in this approximation we can neglect the 
influence of the last term in the Hamiltonian 15]) on the 
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spin-waves of the antiferromagnet and consider the ob- 
tained solutions as excitations of the antiferromagnetic 
part of the heterostructure. 



IV. MAGNETIC ORDERING IN THE METAL 

Now let us consider the influence of the antiferromag- 
netic ordering of the Mott insulator on the magnetic state 
of the metal. In accord with the idea of the spin-wave 
approximation, which considers the spin-wave operators 
&l and &l as describing small deviations from an equilib- 
rium orientation, we can keep the largest terms contained 
in the longitudinal part of the interaction, 



Hu — 2 J SfoSl, 



(30) 



in the last row of Eq. ([8|). This sum contains the zero- 
order terms of the spin- wave approximation S^ = ^e 1 ^ 1 
where Q = (%, ir) [see Eq. ([16]) ]. Thus, we shall consider 
the Hamiltonian 



H„ 



*S 12 fl l+a,U 



laa /,<-! 



*12 12 [4,l x -l,v a ll*<r + Ov^lx-h* 



la 1 X <-1 



J 



12 eiQ1(xa t-i,a a h-h^ 



(31) 



t 



was used. 



where the relation S[_ 1 = | ^ CT aa[ _ 1 a a\^i M 

The parameter we wish to consider is the magnetiza- 
tion in the metal region, 



Mr 



«LT - "Li- 



(32) 



It is convenient to divide the YZ plane into two sublat- 
tices which contain sites with odd and even sums l y + l z 
of the lattice coordinates. The Fourier transformation of 
the electron annihilation operators over the sites of one 
of these sublattices reads 



/ 2 -iql(m) 
\{m) 



m = 1,2, 



where 1(1) and 1(2) label sites of the two sublattices, and 
the wave vector q belongs to the first magnetic Bril- 
louin zone which is half as much as the usual 2D Bril- 
louin zone. These annihilation operators are connected 
with the Fourier-transformed annihilation operators used 
above by the relations 



1 , . 



(33) 



V2 




co/|f| 



FIG. 3: The imaginary parts of Green's functions G^y (ql x l x ) 
(the solid lines) and Gu (ql x l x <y) with a =f (the dashed lines) 
and a =1 (the dash-dotted lines) for q = (0,7r), J = . 23 £ | , 
the temperature T = 0, l x = — 1 (a), and l x — —10 (b). The 
inset in part (b) demonstrates the same curves on an enlarged 
scale. 



where q' is a wave vector in the second magnetic Brillouin 
zone which is connected with the q by the relation q' = 

q+ (±7T,±7r). 

Let us determine the matrix retarded Green's functions 

GW(qtU» = -»0(t) ({a« m iAt),al m%a }) (34) 

with the operator time dependence determined by Hamil- 
tonian (|3"Tj) . The magnetization ([32]) is connected with 
the Fourier transform of this Green's function by the re- 
lation 



M, 



N ^ 
q 



' duj Im [G mm {quil x l x l) - G mm {qtl x l x ])} 



(35) 



where (3 is the inverse temperature and the summation 
over q is over the first magnetic Brillouin zone. 
The equation for Green's function (f3~4")) reads 



i-^ t G{qtl x l' x a) 



G{qtl x l' x cr) 
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FIG. 4: (a) The magnetization M in the metal as a function 
of the distance from the surface l x for J = 0.23|i | and T = 0. 
(b) The absolute value of M in a semi-logarithmic scale for 
T = 0, J = 0-23|t| (filled circles) and J = 0.1|t| (open circles). 



+ t 



G(qt,l x -l,l' x a) 
+G(qi, l x + l,l' x a) 



(36) 



where 



Equation 



T-2 



1 

1 



can be rewritten as 

G(qtlJ' x (r)=G^(qtl x l' x ) 
J 



dt'G^(q,t-t',l x ,-l) 



xG(qt',-l,l' x a), 



(37) 



and finally for the Fourier transform of Green's function 
we find 

G(qcjlJ' x a) = G (0) (q^Q - ^ (0) (qwf a , -1) 
J 



r + -aG(°)(q^,-l,-l)f3 



xr 3 

xG(°)(qw,-l,0 



(38) 



In Eqs. (EB and (25), G<°)(qwU;) is the Green's func- 
tion (|34p in which the operator time dependence and the 



averaging are determined by Hamiltonian ([9]) . It is iden- 
tical to Hamiltonian (|3"Tj) without the last term in the 
right-hand side. Using Eqs. (fT4| and ([33|) we find 

Gil.(q^C) = f — sin(M.)sin(Mi) 
Jo ^ 

I 1 (_i) m + m ' 



(39) 



where q', as before, is the wave vector in the second 
magnetic Brillouin zone, connected with the q by the 
relation q' = q + (±7r, ±7t). 

Equation p8|) looks similar to Eq. (f27|) . However, if in 
the latter equation the multiplier J in the second term of 
the denominator is of the same order of magnitude as the 
energy parameter in Green's function ZM°) (kwOO), in the 
former equation this multiplier is much smaller than the 
energy parameter of G^(qw, —1, —1), the electron band 
halfwidth B, due to the supposition U 3> |t|. There- 
fore surface electronic states do not arise. Nevertheless 
the influence of the second term in the right-hand side 
of Eq. l[55]) is perceptible, especially near the surface, as 
seen in Fig. [3J It is this term which yields the depen- 
dence of Green's function G(ql x l' x a) on the spin projec- 
tion a and leads to a nonzero magnetization (|35j) . As 
seen from Fig. [3Jb) , for wave vectors near the boundary 
of the magnetic Brillouin zone this dependence is percep- 
tible even at large distances from the surface which leads 
to a slow decay of the magnetization with this distance. 
From Eqs. (|3"5|) and (j3"8)) it follows that for a given l x the 
magnetization has opposite signs on the two sublattices 
of the YZ plane. 

The dependence of the magnetization on the distance 
from the surface is shown in Fig. 0J We find that the 
magnetization changes its sign not only on moving from 
one 2D sublattice to the other but also with transfer per- 
pendicular to the surface - the magnetization has antifer- 
romagnetic character. As seen from Fig. 0] this ordering 
penetrates into the metal for distance of tens lattice spac- 
ings. After the initial rapid decrease the magnetization 
decays exponentially with l x with the correlation length 
£ = 5 — 6 lattice spacings [see Fig. Etb)]. The value of £ 
depends only weakly on the ratio J/\t\. 

Let us check the influence of different corrections to 
these results. At first let us consider the contribution of 
the term 6j &io in the spin component Sf in the interac- 
tion (|3"0|) [sec Eq. (fTTj)) ] . The value of this contribution 
can be estimated using the relation 



1 r c 



duj lmDn(kujl x l x ) 



1 



(40) 



Calculations with Green's function obtained in the previ- 
ous section give for zero temperature: (b\ b\o) = 0.0931. 
Thus, the influence of this term reduces to an effective 
decrease of the exchange constant approximately by 20 
percent. Notice that the value of (b\ b\o) is considerably 
grater than the bulk value (&{ ; 6^) = 0.0730. 
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Now let us consider the influence of spin fluctuations in 
the transversal part of the interaction in Hamiltonian (jHJl , 

Hu = JJ2 (s K SLi + > ( 41 ) 

1 

where 

= a i,-i,T a i.-l.lj ^i7-i = a i,-i,j. ai > -1 >T - 
Using the perturbation theory we find 

G'(qnl x l' x a) = G(qnlJ' x a) - G(qnl x , -1 , a) 

xt(qna)G'(qn,~l,l' x a), (42) 

where & (qnl x l' x a) is the Fourier transform of the Mat- 
subara Green's function with the components 

G' mm >(qTl x l x a) = -(Ta qm/x(T (T)a] im , i> ). (43) 

Here T is the chronological operator and the averaging 
and the dependence on the imaginary time r are deter- 
mined with the total Hamiltonian ((SJ. In Eq. (|42]) . the 
integer n stands for the fermion Matsubara frequency 
iuj n = (2n + l)nT and G(qnl x l' x a) is the Green's func- 
tion (|34[) at this frequency. In the Born approximation 
the self-energy in Eq. (|42|) reads 

E mmi (qna) = J 2 — ^ [Atmtsi (qii'OO) 

+ (-l) m+mi IW(q^00)] 

x G mmi (q-q 1 ,n-i/,-l,-l,-CT),(44) 

where D(qvl x l' x ) is the spin-wave Green's function (f24|) 
at the boson Matsubara frequency iu v = 2vkT and fa = 
ma + |(1 - a). 

Notice that Eq. (|42|) has the same structure as Eq. (|37| . 
To estimate the value of the self-energy (j4"4f we can use 
Eq. ([39]) for the electron Green's function and the approx- 
imation for the Green's function D(qi/00) which takes 
into account only the pole of the surface spin wave. Due 
to its dispersion and the integration over k x in Eq. (|39j) 
the Kondo-like divergencies are integrated out. As a re- 
sult it can be seen that the self-energy (|44|) is of the 



order of J 2 / B and is much smaller than the respective 
multiplier J in Eq. (|37p . Thus, the contribution of the 
interaction ([4Tjl can be neglected. 



V. CONCLUSION 

In this paper, we have considered elementary excita- 
tions of the heterostructure of the semi-infinite metal 
and the Mott insulator in the case when the crystals dif- 
fer only in the value of the Hubbard repulsion - it is 
zero in the metal, and it exceeds the critical value for 
the Mott transition in the insulator. At half-filling and 
low temperatures the insulator has the long-range anti- 
ferromagnetic order and its low-lying excitations are spin 
waves. We used the unitary transformation to reduce the 
initial Hamiltonian to a simpler one which describes the 
heterostructure of the metal and the Heisenberg antifer- 
romagnet. At the interface of this heterostructure spins 
of the antiferromagnet interact with spins of electrons in 
the metal. We found that elementary excitations of the 
antiferromagnet are standing spin waves with the disper- 
sion similar to that in the infinite case and the magnon 
mode localized near the surface of the antiferromagnet. 
This mode has the dispersion of the two-dimensional 
spin waves with somewhat increased for a given superex- 
change constant frequency. The description of this mode 
has much in common with the description of localized 
states near a point defect. Analogously to these states 
the mode ejects bulk modes - the standing waves — from 
the region near the surface. The antiferromagnetic order 
of the insulator induces the antiferromagnetic ordering in 
the metal. The magnetization in the metal decreases ex- 
ponentially with distance from the interface. The corre- 
lation length is equal to 5-6 lattice spacings and depends 
only weakly on the parameters of the problem. 
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